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Abstract 

We report on calculations of the formation energies of several [100] and [110] 
oriented step structures on biaxially stressed Si and Ge (001) surfaces. It is shown 
that a novel rebonded [100] oriented single-height step is strongly stabilized by 
compressive strain compared to most well-known step structures. We propose that 
the side walls of "hut" -shaped quantum dots observed in recent experiments on 
SiGe/Si films are made up of these steps. Our calculations provide an explanation 
for the nucleationless growth of shallow mounds, with steps along the [100] and 
[110] directions in low- and high-misfit films, respectively, and for the stability of 
the (105) facets under compressive strain. 
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Strain induced self-assembly is actively being pursued as a technique for the fabrica- 
tion of nanoscale electronic devices and memories that have the potential to bring higher 
speed to information processing and higher areal and volumetric capacity to information 
storage. In the past few years, significant advances have been made on both the techno- 
logical and fundamental aspects of self assembly. On the technological side, it has become 
possible to prepare regular spatial arrays of nanostructures. 

Several basic issues concerning the physical mechanisms involved in the different 
stages of the formation and morphological evolution of the nanostructures are only now 
becoming clear. The initial stages of epitaxial nanostructure formation remains least 
well understood. Of particular interest here, recent experimental work [1|-|4]] has revealed 
that (105) oriented quantum dots in SiGe/Si films grow from shallow precursor mounds 
whose sidewalls are made up of widely spaced steps that are oriented in the [100] or [110] 
directions for low or high Ge compositions, respectively. These observations cannot be 
understood on the basis of a competition between surface energy increase involved in cre- 
ating the walls of the dots and the strain energy reduction through elastic relaxation. Such 
a competition would lead to a nucleation barrier for the formation of the islands, while 
the experiments clearly show that the stepped mounds emerge as a natural instability 
without any such barrier |l].|2[. 

Recently we have demonstrated that nucleationless growth of stepped mounds can 
be understood by including the physics of surface steps, in particular the dependence of 
their formation energy on strain and their interactions ||. In this letter, we calculate 
the parameters that characterize step formation and interactions in order to understand 
the stabilization of [100] or [110] steps on (2x1) reconstructed Si and Ge (001) surfaces 
at different levels of biaxial strain. The strain dependence of the step formation energy 
is determined by the disposition of the atomic bonds in the vicinity of the step-edge; if 
there is significant rebonding such that these bonds are stretched (compressed) relative 
to the bonds in the bulk material, a compressive (tensile) mismatch stress tends to lower 
their formation energy 0. The interactions between steps include the repulsive dipolar 
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interactions J7|, which depend on the local step density, and the non-local attractive 
interactions due to the force monopoles induced by the mismatch strain ||. The surface 
energy of a stepped surface that is oriented at a small angle 9 with respect to the (001) 
direction can be written by including the step formation and dipolar interaction energies 

as 

~ I sin 9 1 ^ 

7 (0, e) = 7o (e) cos 6 + (j3 1 + M\ sin 0| + fc—^jr, (1) 

cos^ 9 

where e is the biaxial surface strain, 70 is the surface energy of the (001) surface, (3\ is the 
formation energy of a step, (3\ is a measure of the sensitivity of the formation energy of a 
step to strain and f3% characterizes the strength of the step-step dipolar interactions. Since 
the attractive monopole induced interactions favor a stepped surface over a flat surface, 
stepped mounds can grow without any nucleation barrier if the surface energy of their 
sidewalls becomes lower than that of the (001) surface ||. It can be seen from (|lj) that 
this indeed happens for stepped surfaces with orientations in the range \9\ < 9 C , where 
$c ~ \J — /3(e) / /?3, if the condition (3{e) = (3\+ Pie < is satisfied. In the present work we 
show by means of atomistic simulations that this condition is satisfied for the mismatch 
strains of interest. 

In the [110] direction, we will only consider the double-height B-type (DB) steps |J; 
these steps are known to be lower in energy than the combination of alternating single- 
height SA and SB steps for vicinal angles that are more than about 1° - 2°. In the [100] 
direction we consider two kinds of steps: (1) zig-zag SA+SB' 100 ' steps proposed by Chadi 
[[]] and experimentally observed by Wu et.al. [PJ and (2) straight steps on which the 
dimer rows of the adjoining terraces arrive at the steps at an angle of 45°. If we focus 
attention on the family of (lOn) surfaces to which the (105) surface belongs, there are 
two distinct low-energy step structures that can be obtained by removing every other 
atom from the step-edge in order to reduce the number of dangling bonds, as shown in 
Fig. |1[ The key difference between the two structures is that in one of them, which we 
call single-height rebonded (SR), there are atoms at the edge of the step bonded to atoms 
on the adjoining terraces, while this type of rebonding is absent in the single-unrebonded 
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(SU) step (refer to Fig. [I]). The structure of the (105) facet made up of SU and SR steps 



correspond to the reconstructions for this surface proposed by Mo et. al. [|Tl|] and Kohr 
and Das Sarma fL2| , respectively. 

We determine the step formation and interaction energies using the following proce- 
dure: the surface energies of various stepped surfaces are first determined as a function of 
applied biaxial strain for both the (lOn) and (lln) families as shown in Fig. |^. The step 
formation and interaction energies are then extracted from the surface energy 7(6, 9) using 
(P. The strain dependent step formation energies for various step structures calculated 



using the Tersoff (T3) |T3[ potential are shown in Fig. |2| for Si and Ge [15|. In all cases, 
the step formation energies can be reasonably approximated by a linear relation assumed 
in ([l]) with small deviations that can be attributed to non-linear effects, particularly at 
larger strains. Calculation of step energies using the Stillinger- Weber (S-W) |14| potential 
also showed almost linear dependence of the formation energies with strain. The param- 
eters (3\ and j3i, obtained by using a linear fit to /3(e), along with the dipolar interaction 
strength /3 3 are given in Table I. 

If we focus attention on the [100] oriented steps in Si, we find that the formation 
energies of the SR steps are lower than those for the SU and SA+SB^ 100 ' steps for com- 
pressive mismatch strains. This can be understood by analyzing the bond lengths at the 
step-edge (refer to Fig. [T]); since there are several stretched bonds in the SR structure, we 
expect its energy to be lower in compressive strains than the SU structure. Although the 
zig-zag SA+SB^ 100 ! steps are energetically lower than the SR steps in the absence of strain, 
Fig. |2| shows that a modest amount of compressive strain (0.3 % in Si) stabilizes the SR 
steps over these steps. Table I also shows that while the formation energies of the SR 
steps on an unstrained surface are larger that the corresponding values for the DB steps, 
the strain sensitivity $1 for the SR steps is about 50% larger in both Si and Ge. This 
implies that the SR steps are preferred over the DB steps at large values of compressive 
strain, with one key difference between Si and Ge. Using both the S-W and T3 potentials, 
we find that the formation energies of the SR steps fall below the energies for DB steps 
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when strains become more than about 1% and about 4% in Si and Ge, respectively (refer 
to the insets in Fig. |2] and Table I). 

The above observations have important implications for the growth of stepped is- 
lands observed in recent experiments []1]-|J. If we assume that the results for Si are 
indicative of the trends in low-misfit films (e ~ l%-2%), it is clear that the [100] SR steps 
should be observed during early stages of island growth since their formation energies are 
lowest among all the cases that we have considered. While we are not aware of a direct 
experimental observation confirming this prediction, our picture is consistent with the 
experiments of Sutter and Lagally [0 who find that in SiGe films with 25% Ge content, 
the stepped mounds are made up of [100] oriented monolayer steps. Experiments that 
can identify the surface structure of these stepped mounds will be invaluable in resolving 
the issue. In distinct contrast to the low-misfit films, experiments on pure Ge || and Si- 
capped Ge [|J islands show that the stepped mounds are made up of [110] oriented steps 

These observations can also be understood on the basis of the formation energy of 



the SR and DB steps in Ge shown in Fig. |2]. Here, the SR steps do not become favorable 
until mismatch strain becomes close to 4% and 5% in the T3 and S-W potentials respec- 
tively; since the strains in these mounds should be less than the mismatch between Si and 
Ge (4.2%), the observation of [110] oriented steps are in agreement with our calculations. 

While a negative step creation energy would lead to formation of surfaces with closely 
spaced steps, repulsive interactions increase as the spacing between the steps become 
smaller. This competition leads to an optimum slope, which can be determined using (JJ) 



as 9* ~ y — /3(e)/3/?3 ||. The existence of such a slope is indeed confirmed by the surface 
energies for Ge shown in Fig. ^|, where at 4% strain the (105) and (115) surfaces are seen 
to be the optimum facets in the (lOn) and (lln) directions respectively. In the case of 
the (lln) surfaces, the (115) surface corresponds to closest spacing of DB steps (the (114) 
and (113) surfaces are reconstructed so that the steps lose their identity). In the presence 
of compressive strains, the (105) surface with SR steps has a smaller surface energy than 
both the high index (lOn) surfaces shown in Fig. ^ and the low index (103) surface. It 
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can also be seen from Fig. ^| that the surface energy of the (105) surface in Ge lies below 
the energy of the (115) surface throughout the range of interest, while it falls below the 
energy of the (117) at a strain of about 3%. In the case of Si, we find a similar trend 
with the (105) surface falling below the (117) surface close to 1% . These observations 
indicate that the (105) facet made up of SR steps is stabilized by compressive strains in 
both low- and high-misfit SiGe films. On the other hand, the (105) surface with the SU 
steps proposed by Mo et.al. is stabilized by tensile rather than compressive stresses since 
the bonds at the step-edge are in compression (refer to Fig. [I]). While the results in the 
present work are obtained using empirical potentials , we have verified by means of ab 
initio simulations that the atoms at the edges of the SR steps in the (105) surface are 
indeed in a stretched state. The details of these calculations will be published elsewhere. 

In summary, the nucleationless growth of stepped mounds can be understood on the 
basis of strain dependence of [100] and [110] oriented steps on (001) surfaces. The rebond- 
ing at the edge plays an important role in stabilizing stepped surfaces in the presence of 
strain. The competition between repulsive step interactions and strain induced lowering 
of step formation energies leads to optimum low-energy orientations such as (105). Fur- 
ther work that analyzes the effect of surface segregation on step formation energies will be 
invaluable in gaining deeper insight into early stages of quantum dot formation in SiGe 
and other alloys. 

Acknowledgments : The research support of the National Science Foundation through grant 
CMS-0093714 and the Brown University MRSEC Program, under award DMR-0079964, is grate- 
fully acknowledged. 



6 



References 

1. P. Sutter and M. G. Lagally, Phys. Rev. Lett. 84, 4637 (2000). 

2. R. M. Tromp, F. M. Ross and M. C. Reuter, Phys. Rev. Lett. 84, 4641 (2000). 

3. A. Vailionis, B. Cho, G. Glass, P. Desjardins, D. G. Cahill and J. E. Greene, Phys. 
Rev. Lett. 85, 3672 (2000). 

4. A. Rastelli, M. Kummer and H. von Kanel, Phys. Rev. Lett. 87, 6101 (2001). 

5. V. B. Shenoy and L. B. Freund, A continuum description of the energetics and evolu- 
tion of stepped surfaces in strained nanostructures, J. Mech. Phys. Sols, (in press); 
A copy of this article can be found at \http:// 'arXiv.org/abs/cond-mat/02035 14\ 

6. Y. H. Xie, G. H. Gilmer, C. Roland, P. J. Silverman, S. K. Buratto, J. Y. Cheng, 
E. A. Fitzgerald, A. R. Kortan, S. Schuppler, M. A. Marcus and P. H. Citrin, 73 
3006 (1994); C. Roland, MRS Bulletin 21, 27 (1996). 

7. V. I. Marchenko and A. Ya. Parshin, Sov. Phys. JETP 52, 129 (1980). 

8. J. Tersoff, Y. H. Phang, Z. Y. Zhang and M. G. Lagally, Phys. Rev. Lett. 75, 2730 
(1995). 

9. D. J. Chadi, Phys. Rev. Lett. 59, 1691 (1987). 

10. F. Wu, S. G. Jaloviar, D. E. Savage and M. G. Lagally, Phys. Rev. Lett. 71, 4190 
(1993). 

11. Y. W. Mo, D. E. Savage, B. S. Swartzentruber and M. G. Lagally, Phys. Rev. Lett. 
65, 1020 (1990). 

12. K. E. Khor and S. Das Sarma, J. Vac. Sci. Technol. B, 15 1051 (1997). 

13. J. Tersoff, Phys. Rev. B, 38, 9902 (1988); Phys. Rev. B, 39 5566 (1989). 



7 



14. F. H. Stillinger and T. A. Weber, Phys. Rev. B, 31, 5262 (1985); K. Ding and 
H. C. Andersen, Phys. Rev. B, 34 6987 (1986). 

15. In the case of SR steps there should be an additional term in ([!]) that accounts 
for interaction of the monopoles due to the discontinuity of surface stress at the 
step-edge. We have ignored this term in our analysis of vicinal surfaces with 9 > 3° 
since the magnitude of this term becomes significant for small vicinal angles (6 < 

1°)- 

16. From the STM images of the stepped mounds in || we have inferred that the [110] 
oriented steps are of the double-height type. 



8 



Table I: Step formation energies Pi, their derivatives with respect to strain Pi, and the 
step-step dipolar interaction coefficient /3 3 for Si and Ge (with values for Ge enclosed in 
parentheses) calculated using the Stillinger- Weber (S-W) and Tersoff (T3) potentials. In 
order to allow a direct comparison of the formation energies of different step structures, 
Pi is measured relative to the formation energy of the DB step calculated using the T3 
potential. All quantities are expressed in meV/ A 2 . 





SR 




DB 




SU 


SA+SBl 100 ! 




S-W 


T3 


S-W 


T3 


T3 


T3 


Pi 


19.9(29.9) 


5.8(10.1) 


12.1(14.8) 


0.0(0.0) 


46.3(38.6) 


5.98(0.96) 


Pi 


741(738) 


1046(781) 


440(459) 


667(501) 


-93(-89) 


523(385) 


& 


169(186) 


217(143) 


239(189) 


307(229) 


2(11) 


585(382) 



9 



Figure Captions 



Figure 1: The structure of the [100] oriented single-height rebonded (SR) step (top) and 
the unrebonded (SU) step (bottom) on a (2x1) reconstructed (001) surface. The dashed 
line denotes the edge of the step, where the dimer orientation undergoes a 90° rotation. 
Atoms are colored according to the number of dangling bonds(b) before surface and step- 
edge reconstructions: red=2b, green=lb, blue=0b. The numbers represent the stretching 
(in % ) of some of the bonds relative to the bond-length in the bulk for the Si(109) surface 
relaxed using the T3 potential. Note that the rebonding of the atoms across the step edge 
leads to stretching of several bonds in the SR step. 



Figure 2: Surface energy of stepped surfaces (in meV/ A 2 ) consisting of [100]SR(red) and 
[110]DB(black) steps as a function of the biaxial strain e for Ge, computed with the T3 
potential. The vicinal angles and surface orientations are indicated in the figure. The 
insets show the strain dependent formation energy /3(e) (in meV/A 2 ) for the three types 
of steps, DB (black squares), SR(red circles) and SA+SBt 100 ' (blue triangles) in Si and Ge. 
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Figure 1: Shenoy, Ciobanu and Freund 
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Figure 2: Shenoy, Ciobanu and Freund 
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